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Abstract 

The Laser Interferometer Gravitational Wave Observatory (LIGO) is a net- 
work of three detectors built to detect local perturbations in the space-time 
metric from astrophysical sources. These detectors, two in Hanford, WA and 
one in Livingston, LA, are power-recycled Fabry-Perot Michelson interferom- 
eters. In their fifth science run (S5), between November 2005 and October 
2007, these detectors accumulated one year of triple coincident data while 
operating at their designed sensitivity. In this paper, we describe the calibra- 
tion of the instruments in the S5 data set, including measurement techniques 
and uncertainty estimation. 

Keywords: Interferometer, Calibration, Control Systems, Gravitational 
Waves 

PACS: 07.60.Ly, 07.05.Dz, 04.30.-w, 04.80.Nn 



1. Introduction 

The Laser Interferometer Gravitational Wave Observatory (LIGO) is a 
network of three detectors built in the United States to detect local pertur- 
bations in the space-time metric from astrophysical sources. These distant 
sources, including binary black hole or neutron star coalescences, asymmet- 
ric rapidly spinning neutron stars, and supernovae are expected to produce 
time-dependent strain h{t) observable by the interferometer array [37i [21j. 



6 



The detectors, two in Hanford, WA (HI and H2) and one in Livingston, 
LA (LI), are power-recycled Fabry-Perot Michelson interferometers. The 
optical layout of the interferometers is shown in Figure [T] The perpendicular 
Fabry-Perot arm cavities of the Michelson, each of length L = 3995 m for 
HI and LI (L = 2009 m for H2), are composed of 10 kg optics or "test 
masses" suspended as pendula. Light reflected from the input port of the 
Michelson is recycled with an additional suspended optic forming a power 
recycling cavity. Each interferometer uses a Nd:YAG laser (A = 1064 nm, 
or / = 282 THz), whose phase is modulated at several frequencies such 
that a Pound-Drever-Hall style control scheme [131 EI] can be used to hold 
the arm cavities and power recycling cavity in resonance. Figure [2] shows 
a schematic of the suspension system for a given optic and electro-magnetic 
coil-actuators (paired with magnets secured on the rear face of the optic) used 
to control its motion. Further details of the interferometer conflguration are 
described in p]. 




Figure 1: Schematic optical layout of the LIGO interferometers. 

During the flfth LIGO science run (S5), these detectors accumulated ap- 
proximately one year (368.84 days) of triple coincidence data near their de- 
signed sensitivity between Nov 4 2005 and Oct 1 2007 (GPS time 815097613 
through 875232014). The best sensitivity (strain amplitude spectral density) 
for each detector and an example sensitivity curve used to guide the design 
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Figure 2: A schematic of the LIGO optic suspensions for S5. The actuation force is 
provided by the coil actuators (mounted to the support structure) which act upon the 
magnets secured directly on the rear face of the optic. 



for the 4 km detectors ^ are shown in Figure [3| As a figure of merit of 
the sensitivity over time, we integrate the power spectral density using a 
matched- filter template describing a binary neutron star (1.4-1.4 solar mass) 
coalescence over which angle and orientation have been averaged. This met- 
ric produces a predicted range out to which we may see such a source with 
signal-to-noise ratio of 8 (see [4J for details) . Figure [i] illustrates the daily 
median of this range over the course of the science run. 

Differential displacement of the interferometer's end test masses is mea- 
sured by precisely monitoring the differential phase between light returned by 
each Fabry-Perot arm cavity using a Pound-Drever-Hall error signal. When 
the interferometer is under servo control, this error signal euif) is propor- 
tional to a differential arm (DARM) length change, AL^xtif) caused by the 
end test mass displacement such that, 

ALUf) = Rdf) enif) (1) 

where the change in length AL^xt is the sum of the interferometer's response 
to the astrophysical signal and other differential noise sources. 

The quantity Rhif) is a complex function in the frequency-domain known 
as the "length response function." In this paper, we provide a complete de- 
scription of a frequency-domain model of the length response function used 
for each detector in the S5 data set. Table [T] summarizes the uncertainty in 
our model of i?L(/), broken up into magnitude and phase of the complex 
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Best Displacement Sensitivity of LIGO Interferometers 




Frequency (Hz) 

Figure 3: The best displacement sensitivity, expressed as equivalent displacement noise, 
for each interferometer during the S5 science run, and expected total noise in LIGO's first 
4 km interferometers. 

function, and separated into three frequency bands. Each value is the esti- 
mated 68% confidence interval (one sigma) across the band for the entire 2 
calendar-year science run. 



Table 1: Summary of band-limited response function errors for the S5 science run. 





RLif) Mag 
40-2000 Hz 


nitude Ei 

2-4 kHz 


[■ror (%) 

4-6 kHz 


Riif) Phj 
40-2000 Hz 


ise Error 

2-4 kHz 


(Deg) 

4-6 kHz 


HI 
H2 
LI 


10.4 
10.1 
14.4 


15.4 
11.2 
13.9 


24.2 
16.3 
13.8 


4.5 
3.0 
4.2 


4.9 
1.8 
3.6 


5.8 
2.0 
3.3 



In Section [2| we describe the model used for all LIGO interferometers 
which divides a given interferometer into three major subsystems - sens- 
ing, digital control, and actuation - and includes a detailed description of 
the important components of each subsystem. Measurements of these com- 
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Figure 4: Daily median of detector sensitivity during S5 to a 1.4-1.4 solar mass compact 
binary system averaged over angle and orientation. Dashed lines indicate the times during 
which the representative spectra in Figure [3] were taken. Large variations in detector 
sensitivity are due to upgrades or hardware problems. 



ponents along with corresponding uncertainties are presented in Section |3| 
Finally, the response function, i?L(/), is developed from the subsystems and 
the uncertainty in each subsystem are combined in Section [4] to form the 
total uncertainty estimate as seen in Table [T] 

Gravitational wave data analysis is performed on a signal proportional 
to strain generated in the time domain from eD{t) and a convolution kernel, 

RL{t-t'), 

h{t) = ^J RL{t-t') eoit') dt', (2) 

developed from the parameters of the length response function. The produc- 
tion of the time-domain convolution kernel, RL{t — t')^ from the frequency- 
domain model, RlH)^ and the associated additional uncertainty are dis- 
cussed in detail in j33l [22] . 
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2. Model 



Astrophysical gravitational wave strain h{f) detected by the interferom- 
eters contains source information including wave forms h^^x{f)^ azimuthal 
angle 0, polar angle 9^ and orientation (or polarization angle) ip (see Figure 
[5]). The amplitude of the wave's projection into the interferometer basis is 
described by 



where Fx,+ are the antenna response of the detectors and /ix,+ are the wave 
amplitudes in the "cross" and "plus" polarizations of the local metric per- 
turbations hj^jy in the transverse-traceless gauge [371 El E] • 



Figure 5: A schematic of the coordinates used both in the interferometer basis, and in the 
incoming plane wave basis. The Euler angles 6>, and i/j are as defined in PI E], except 
here they are shown relative to the detector frame rather than the equatorial frame. 

We model each interferometer's repsonse to an optimally-oriented {9 = 
(/) = ip = 0)^ plus-polarized wave form using the long wavelength approxima- 
tion. The approximation is valid between 40 and 6000 Hz, and has associated 
uncertainty of at most 2% [29]. From this reference model, the detector re- 
sponse to an arbitrary waveform, orientation, and polarization angle may be 
calculated analytically [371 [9l [8]. I^i the long wavelength approximation, the 
strain amplitude, h{f)^ in Fabry- Perot arm lengths of the interferometer is 



h{f) = Fx {e, 0, ^)/ix (/) + F+{e, (p, ij)h+{f) 



(3) 




hif) 



'ext 



(4) 



L L 
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Feedback control systems are used to hold the interferometer in a regime 
where the digital error signal, eD(/), is linearly related to the DARM length, 
ALext^ (as in Eq. [T]) and hence to the gravitational wave strain, h{f). We 
model this control loop as a single-input, single-output control loop depicted 
in Figure [6| 



AL 



ext' 



y(t) C 



L 



bo' H 
^ I Q 

^ I H 

^ ' > 
! 



A 





— ► 

r 


D 



D 



Figure 6: The single-input, single-output model of the control loop of differential end test 
mass motion. The interferometer senses and digitizes a change in DARM length, AL^xt 
according to ^{t)CL{f)^ the result of which is the digital error signal which is then fed 
back through a set of digital filters D{f)^ and converted to analog control via the actuation 
function of the end test masses A{f) . 

The loop contains three major subsystems. First is the length sensing 
function, CL{f^t)^ which describes how the interferometer responds to dif- 
ferential changes in arm lengths and how that response is digitized. This 
function is separated into a frequency-dependent function which may 

have some slow time dependence captured by a factor j{t). D{f) is a set 
of digital filters, used to shape the loop error signal into a control signal. 
The remaining subsystem is the actuation function, ^(/), which describes 
how the test masses physically respond to the digital control signal. We as- 
sume linear relationships between all subsystems, such that any subsystem 
(and internal components) may be defined by the ratio of output over input 
signals. 

The product of frequency-dependent subsystems inside the control loop 
is the "open loop transfer function" 

GLif) = CL{f)D{f)A{f). (5) 

Using the above model, we derive the length response function, RL{f,t), in 
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terms of these functions to be 



The remainder of this section describes the components of each subsystem 
in the control loop. 

2.1. Sensing Function 

The length sensing function, CL{f^t)^ describes the transfer function be- 
tween the residual change in DARM length, AL(/), and the digital error 
signal, eD(/), 

Cdf^t) = 7^1^. (7) 

It is important to note that this linear relationship between the DARM length 
change and the digital error signal only applies when the detector is under 
control of the feedback loop: in Eq. [Tj AL(/) is the residual external DARM 
length change, ALextif)^ after the controlled length change, AL^(/), is ap- 
plied. The sensing function has several components (shown in Figure [?]) 
which are treated independently, 

CL{f,t) = 7{t) X JCc X [CFp{f) X ADCif)]. (8) 

The constant, /Cc, which holds all frequency-independent scaling factors, has 
dimensions of digital counts of error signal per unit change in DARM length. 
The remaining terms in Eq. [8] are dimensionless, including time dependence, 
treated independently in the coefficient j{t). 

The change in each arm cavity length L affects the phase of the laser's 
electric field returning from the cavity. On resonance, the transfer function 
between the change in electric field phase refiected by the cavity input mirror 
$(/) and a change in cavity length is 

^ 27rl re(l-r|)sin(27r/L/c) e-'^^f^/^ 
^"'^^^ Ar, (l-r^re) 27r/L/c I - nree-^^'f^l'^ ^ ' 

where A is the laser wavelength, Tc = (re — Tj)/(1 — riVe) is the on-resonance 
Fabry-Perot arm cavity reflectivity, rj and are the amplitude reflectivity 
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Figure 7: Schematic breakdown of the sensing function CL{f,t). Internal to the LIGO 
Scientific Collaboration, the digital signal ejj is often colloquially referred to by its digital 
"channel" name DARM_ERR. From left to right, CfpU) ^ ^fp + ^fp is the arm 
cavity transfer function; a{t) is the time-dependent variation of the interferometer's input 
laser power and optical gain; JCc is the scaling coefficient which absorbs all constants 
including the input laser power, optical gain, the quantum efficiency of the photodiodes, 
the impedance of the photodiode circuitry, and the analog-to-digital gain; ADC{f) is the 
frequency dependence of the analog to digital conversion; and P{t) is the digital factor 
which compensates for the analog change a{t). The compensation is not perfect, therefore 
the factor ^{t) = a{t)l3{t) represents the residual variation. 



of the input and end test masses, and c is the speed of light. In the fre- 
quency band considered for analysis, where 40 Hz < / < 6 kHz <C 2c/L, the 
frequency-dependence of Hpp^f) is approximated by a simple "cavity pole" 
transfer function, 

_ H^pU « C/2L) ^ 1 

""''^^^ = HMO) ITTf ' 

where = c (1 - r,re)/47rLy^ [ig [29l EQl EH [35] . 

The LIGO detectors use a Pound-Drever-Hall detection scheme to ex- 
tract this phase information from the arm cavities, which is recombined at 
the beam splitter. The laser electric field input into the interferometer is 
phase-modulated at uJm/'^TT =25 MHz, which effectively splits the field into 
a "carrier" field with the original laser frequency, f], and upper and lower 
"sideband" fields with frequency f^zbcj^. The sideband fields resonate in the 
power recycling cavity but are anti-resonant in the arm cavities, and there- 
fore, unlike the carrier field, experience no phase change from the arm cavity 
length variation. The Michelson is set up with a fixed asymmetry such that, 
at the anti-symmetric port, the carrier field is held on a dark fringe and the 
sideband fields are not. In this setup, when the arm cavity lengths change 
differentially, the carrier field moves away from the dark fringe, mixes with 
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the sideband field at the antisymmetric port, and a beat signal at uj^ is 
generated. 

The power of the mixed field at the antisymmetric port (in Watts) is 
sensed by four photodiodes. The photocurrent from these diodes is con- 
verted to voltage, and then demodulated at 25 MHz. This voltage signal 
(and therefore the change in DARM length) is proportional to power of the 
input laser field, the "optical gain" (the product of Bessel functions of mod- 
ulation strength, the recycling cavity gain, the transmission of the sidebands 
into the antisymmetric port from the Michelson asymmetry, the refiectivity 
of the arm cavities for the carrier), the quantum efl&ciency of the photodi- 
odes, and the impedance of the photodiode circuitry [3ll [IHl |34l [35]. The 
demodulated voltage from the photodiodes is whitened, and anti-aliased with 
analog circuitry and then digitized by an analog-to-digital converter which 
scales the voltage to digital counts. The frequency dependence of the anti- 
aliasing filters and digitization process is folded into the function ADC{f). 
We absorb all proportionality and dimensions of this process into the single 
constant, /Cc, having dimensions of digital counts per meter of DARM test 
mass motion. 

The optical gain is time-dependent because small, low-frequency (/ <C 
40 Hz) alignment and thermal lensing fiuctuations in the resonant cavities 
change the carrier and sideband field amplitudes. The input laser power 
may also fiuctuate from similar alignment and thermal effects. We represent 
these variations with a coefficient, a{t). The input power, along with the 
carrier and sideband power stored in the cavities, are monitored by several 
independent photodiodes. Their signals are also digitized and combined to 
form a coefficient, /5(t), used to digitally compensate for the time-dependent 
variations. The compensated anti-symmetric port signal forms the error 
signal for the DARM control loop, eoif)- The sensing function therefore 
depends on both time and frequency, but can be separated into independent 
components Cl(/, t) = ^{t)CL{f)^ where 

CLif) Cpp{f)xADC{f) = [HUf) + mpif)]xADC{f) (11) 

and 

7(t) = a{t)/3{t) (12) 

is the scale factor of order unity accounting for the residual time dependence 
after compensation. 
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2.2. Digital Filters 

The digital filters, D{f)^ are known functions in the model. These filters 
are used to shape the digital DARM control loop error signal, euif) (in digital 
counts proportional to displacement) into a digital control signal, s^if) (in 
digital counts proportional to force), 

m = j47t (13) 

eD[f) 

Over the course of the science run, discrete changes are made to the digital 
filters, D{f)^ to improve the performance and stability of the detector (four 
times in the Hanford interferometers, three in Livingston). These changes 
significantly alter the frequency-dependence of the DARM control loop, and 
hence aflFect the overall response function of the interferometer. We divide 
the run into "epochs" defined by these changes. 

Note that the digital filter component does not include all digital filters 
in the DARM loop. Both the sensing function and the actuation function 
contain digital filters, but their frequency dependence is either negligible in 
the measurement band, only important in a very narrow frequency range, or 
are compensating for analog circuitry whose product with the digital filters 
form a unity transfer function. We include these filters in their respective 
sub-systems for completeness. 

2.3. Actuation function 

The actuation function A{f) is defined by the transfer function between 
the digital control signal, s^if)^ and the physical motion imposed on the end 
test masses by the control loop, AL^, 

A{f) = (14) 

and has units of end test mass displacement in meters per count of digital 
control signal. We describe the actuation function as a linear combination of 
functions for each test mass, 

Aif) = eA%f) + eAyif). (15) 

where (^^'^ are known digital coefficients of order unity, roughly equivalent, 
but opposite in sign. Once split, the control signal ffows through each com- 
ponent to the end test masses as shown in Figure [Sj 
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Figure 8: Schematic breakdown of the signal flow through the actuation function for the 
X arm A^{f). The digital signal soif) is colloquially referred to by its "channel" name 
DARM.CTRL. From right to left, is the fraction of the digital control signal sent to the 
X arm; D\{f) are digital fllters; DAC{f) is the frequency dependence of digital to analog 
conversion; Kj^ is the scaling coeflicient proportional to the digital-to-analog gain, the 
gain of resistance circuitry which converts voltage to current, the gain of the coil actuators 
which convert current to magnetic force, and the force-to-displacement transfer function 
gain; and is the frequency dependence of the force-to-displacement transfer function. 



For each arm, the digitally split control signal passes through digital 
suspension filters, and is converted from digital signal to an analog 

voltage via the digital to analog conversion element, DAC{f)^ which includes 
analog anti-imaging circuitry. The resulting voltage passes through a resis- 
tance circuit converting it into current, and is sent to the coil actuators which 
convert the current into force on the magnets attached to the end test mass. 
The suspended test mass is displaced according to the force-to-displacement 
transfer function, changing each arm cavity length, AL^(/). The arm's 

scaling coefficients, /C^, absorb all dimensions and frequency-independent 
factors in the actuation path. This includes the digital-to-analog gain, the 
gain of the resistance circuitry, the gain of the coil actuators, and the force- 
to-displacement transfer function scale factor. In summary, we express the 
individual end test mass actuation functions in Eq. [15] as 

A^n/) = K:Yx[DY{f)xDAC{f)xP^'y{f)\. (le) 

The actuation coefficients, /C^, scale the arbitrary counts of digital excitation 
force into meters of test mass motion. The remaining terms in Eq. [16] are 
dimensionless. 

The suspended test mass can be treated as a pendulum driven by the coil 
actuators (see Figure [2]). The force-to-displacement transfer function for the 
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center of mass of a pendulum, Pcmif), is 



Pcmif) oc 



[jri 



f-P 



(17) 



Qcm 



where /q^ and Q^'^ are the frequency and quahty factor of the pendulum. A 
rigid body resonant mode akin to the fundamental mode of a cylindrical plate 
[27] (see Figure |9]) known as the "drumhead" mode is also included in the 
force-to-displacement model. Its radially-symmetric shape, excited by the 
actuators, lies directly in the optical path and amplifies the cavity's response 
to the length control signal above a several kHz [T2j. We approximate the 
effects of the resonance by multiplying Pcm{f) by an additional pendulum 
transfer function, defined by frequency, /q^, and quality factor, Q^^. 
The total force-to-displacement transfer function is 



Figure 9: Physical shape of the end test mass drumhead internal resonance. Left: Cartoon, 
edge-on view of the fundamental mode of a cylindrical plate [27J. Right: Three dimensional 
modal shape of the drumhead resonance from finite element analysis of a cylinder with 
dimensions similar to the LIGO test masses [12j. 

The digital suspension filters, are between the split control signal 

and the digital-to-analog converter. Their purpose is to remove control signal 
in narrow frequency ranges around the frequencies of other in-band, non- 
axisymmetric, rigid-body resonant modes of the test masses that are excited 
by the actuation forces [12j, and to reduce the coupling between DARM 
length motion and angular motion of the test mass. 



(18) 
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3. Measurements 



Each subsystem of the response function Rhif) is developed using mea- 
surements of key parameters in their modeled frequency dependence and 
their scaling coefficients. The digital filter subsystem is completely known; 
its frequency dependence and scaling coefficient are simply folded into the 
model of the response function. The parameters of the frequency-dependent 
portions of the sensing and actuation subsystems may be obtained precisely 
by direct measurement or are known from digital quantities and/or design 
schematics. As such, these parameters' measurements will only be briefly 
discussed. 

The detector's sensing function behaves in a non-linear fashion when un- 
controlled, therefore we may only infer the linear model's scaling coeflicient, 
Kc-) from measurements of the detectors under closed control loops. We 
infer that the remaining magnitude ratio between our model and measure- 
ments of the open loop transfer function as the sensing coefficient 
Kc (where fuGF is the unity gain frequency of the DARM control loop). 
Other than the known frequency-independent magnitude of D{f)^ the open 
loop gain model's magnitude is set by the actuation scale factor, JCa- This 
makes it a crucial measurement in our model because it sets the frequency- 
independent magnitude of the entire response function. Measurements of the 
open loop transfer function over the entire gravitational wave frequency band 
are used to confirm that we have modeled the correct frequency dependence 
of all subsystems. Finally, measurements of j{t) track the time dependence 
of the response function. The details of these measurements and respective 
uncertainty estimates are described below. 

3.1. Actuation Function 

The components of each arm's actuation function, /C^'^, DAC{f) and 
P^^y(^f) are measured independently in a given detector. As with D{f)^ 
both (^^'^ and D^^(f) are digital functions included in the model without 
uncertainty. 

3.1.1. Actuation Scaling Coefficients, 

The standard method for determining the actuation coefficients, /CJ^^ 
used for the fifth science run is an interferometric method known as the 
''free-swinging Michelson" technique; a culmination of several measurements 
with the interferometer in non-standard configurations. The method uses 
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the interferometer's well-known Nd:YaG laser wavelength (A = 1064.1 ±0.1 
nm, [241136] ) as the calibrated length reference while using the test mass' coil 
actuators to cause a length change. Details of the technique are described in 
[Appendix A[ The actuation coefficient is measured using this method many 
times for each optic in each interferometer over the course of the science 
run, and their mean used as the actuation scaling coefficient for all model 
epochs. Table [2] summarizes the actuation coefficients, /C^'^, for the three 
interferometers in the fifth science run, using free-swinging techniques. 



Table 2: Summary of the actuation scaling coefficients measured during S5. These single 
numbers are formed by the mean of each measurement's median {JC^^)j (6 for each end 
test mass in HI, 5 in H2, and 14 and 15 for the X and Y test masses, respectively in LI). 
Only statistical uncertainty is reported here; systematic uncertainty is folded the total 
uncertainty of the actuation function. 





/C^ (nm/ct) 


/C^ (nm/ct) 


HI 
H2 
LI 


0.847 ±0.024 
0.934 ± 0.022 
0.433 ± 0.039 


0.871 ±0.019 
0.958 ± 0.034 
0.415 ±0.034 



3.1.2. Force-to- Displacement Transfer Function^ P{f) 

Each test mass coil actuator system is equipped with an optical position 
sensor system that consists of an infrared LED emitter aimed at a small 
photodiode mounted in the coil actuator, and a mechanical "fiag" attached 
to the magnet on the optic that cuts through the beam. From amplitude 
spectral densities of these sensor signals while the optic is free-swinging, the 
frequency of each center-of-mass transfer function, /q^ is measured with 
negligible uncertainty. The quality factors, Q^^, depend on the amount of 
local damping applied to suspension, but are estimated from driven transfer 
functions. The uncertainty of this estimation, though large, has little effect 
on the center-of-mass transfer function in the frequency band of interest 
and is ignored. Table [3] shows the results for the center of mass force-to- 
displacement transfer function. The drumhead frequency, /q^ for each test 
mass in the Hanford and Livingston detectors have been measured to be 9.20 
kHz and 9.26 kHz respectively with Q^^ ^ 10^ [121 EH], where again though 
the uncertainty in these parameters may be large, it has little effect in band 
and is ignored. 
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Table 3: Summary of pendulum frequencies, /q^, and quality factors, Q^^, used to com- 
pose models of each interferometer's center-of-mass pendulum transfer functions in S5. 





X End Test Mass 


Y End Test Mass 




/o=- (Hz) 




/o^-(Hz) 




HI 


0.767 


10 


0.761 


10 


H2 


0.749 


10 


0.764 


10 


LI 


0.766 


100 


0.756 


100 



3.1.3. Digital to Analog Conversion^ DAC{f) 

The digital to analog conversion model DAC{f)^ includes the effects of 
the finite sample-and-hold method used to convert digital signal to an analog 
voltage, the analog anti-imaging filter, measured residual frequency depen- 
dence from imperfect digital compensation of analog de-whitening, and the 
time delay arising from computation and signal travel time. 

We use the standard model for the sample and hold of the digital to 
analog converter ||28l [TT] 

H,{f) = sinc[(27r/)/(2/,)] e-^(2-/)/(2/3), (19) 

where the sample frequency fg = 16384 Hz is used in all detectors. 

The same analog anti-image filter is used for each of the four coils on the 
test mass. They are analog, third-order, Chebyshev low-pass filters with 0.5 
dB passband ripple whose corner frequency is at 7.5 kHz and 8.1 kHz for the 
Hanford and Livingston detectors, respectively, and modeled as such in the 
DAC{f) transfer function. We also include residuals measured between the 
modeled anti-imaging filter and its analog counterpart. 

For a given end test mass, there is a complementary pair of digital and 
analog whitening filters for each of the four coil actuators. A comparison 
between the digital compensation and the real analog electronics has shown 
non-negligible, frequency-dependent residuals. We measure the residuals for 
all four coils in each test mass by taking the ratio of transfer functions be- 
tween a digital excitation and the analog output of the whitening filters with 
the digital filters on and off. We include the average residual of the four coils 
in our model. 

A detailed analysis of the digital time delay in the digital-to-analog con- 
version has been performed elsewhere [lOj. For the actuation model we es- 
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timate the time delay from our model of the open loop transfer function 
(attributing all residual delay in the loop to the actuation function), and 
assign a fixed delay to each epoch. 

3.1.4' Actuation Uncertainty^ a a 

The digital suspension filters, DA{f)-> have well-known digital transfer 
functions, which are included in the model without an uncertainty. The 
model of force-to-displacement transfer function, and digital-to-analog 

conversion, DAC{f)^ are derived from quantities with negligible uncertainty. 
Hence, the uncertainty estimate for the actuation function is derived entirely 
from measurements of the actuation scaling coefficient, Ka- 

The actuation coefficient is measured using a series of complex transfer 
functions taken to be frequency independent as described in [Appendix A 



We take advantage of this fact by estimating the frequency-independent un- 
certainty in the overall actuation function from the statistical uncertainty of 
all free-swinging Michelson measurements. For magnitude, we include a sys- 
tematic uncertainty originating from an incomplete model of the actuation 
frequency-dependence, such that the total actuation uncertainty is 



< = ^L.- (21) 



\A\J \\1Ca\J \{r/a) 

2 _ 2 

The statistical uncertainties, cr\]c^\/\lCA\ and cr^^^ , are the quadrature sum 
of the scaling coefficient uncertainty from each test mass, as measured by the 
free-swinging Michelson technique. For each optic's coefficient, we estimate 
the uncertainty by taking the larger value of either the standard deviation 
of all measurement medians, or the mean of all measurement uncertainties 
divided by the square root of the number of frequency points in a given mea- 
surement. These two numbers should be roughly the same if the measured 
quantity followed a Gaussian distribution around some real mean value and 
stationary in time. For all optics, in all interferometers, in both magnitude 
and phase, these two quantities are not similar, implying that the measure- 
ments do not arise from a parent Gaussian distribution. We attribute this to 
the quantity changing over time, or a systematic error in our measurement 
technique that varies with time. Later studies of the free-swinging Michelson 
technique have revealed that the probable source of this time variation is our 
assumption that the optical gain of the simple Michelson remains constant 



over the measurement suite (see Appendix A). 
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We have folded in an additional a{^ria)/ {r / a) = 4% systematic error in 
magnitude for the Hanford detectors only. This correction results from the 
following systematic difference between the Hanford and Livingston free- 
swinging Michelson measurement setup. Analog suspension filters, common 
to all detectors, are used to increase the dynamic range of the coil actuators 
during initial control of the test masses. When optic motions are sufficiently 
small enough to keep the cavity arms on resonance, they are turned off and 
left off as the detectors approach designed sensitivity [2114]. These additional 
suspension filters were left in place for the Hanford measurements in order to 
obtain better signal-to-noise ratios for the driven transfer functions described 
in [Appendix A[ The filters' color had been compensated with digital filters, 
but the average residual frequency dependence is roughly 4% for both end 
test masses in HI and H2. 

The total uncertainty for each interferometer's actuation function, as de- 



scribed in Eq. [21} is shown in Figure [10) These estimates include statistical 
and known systematic uncertainties. To investigate potential unknown sys- 
tematic uncertainties in the actuation functions we applied two fundamen- 
tally different calibration methods. The results of these investigations are 
described in section 5. 



i 0.14 
■ 0.12 

; 0.1 

1 0.08 
■0.06 
: 0.04 
: 0.02 
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10-^ 



Figure 10: Summary of the actuation uncertainty for all detectors in S5. 
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3.2. Sensing Function 

The components of the sensing function, /Cc, Cppif)^ and ADC{f) are 



described in §2.1[ The frequency-dependent components are developed from 
measured parameters with neghgible uncertainties, and Kc is obtained as 
described above. The techniques used to obtain the parameters are described 
below. 

3.2.1. Sensing Scaling Coefficient^ Kc 

In principle, the scaling coefficient Kc is also composed of many inde- 
pendently measurable parameters as described in §2.1[ In practice, these 
components (specifically components of the optical gain) are difficult to mea- 
sure independently as the interferometer must be controlled into the linear 
regime before precise measurements can be made. The scaling coefficient for 
the other subsystems are either measured (in the actuation) or known (in 
the digital filters). We take advantage of this by developing the remainder 
of sensing subsystem (i.e. its frequency-dependence), forming the frequency- 
dependent loop model scaled by the measured actuation and known digital 
filter gain, and assume the remaining gain difference between a measurement 
of open loop transfer function and the model is entirely the sensing scale 
factor. Results will be discussed in §3.3[ 

3.2.2. Fabry-Perot Cavity Response^ Cppif) 

Our model of the Fabry-Perot Michelson frequency response is the sum of 



the response from each arm as in Eq. [11} Using the single pole approximation 
(Eq. 10), the frequency response of each arm cavity Hgp{f) can be calculated 



explicitly using a single measured quantity, the cavity pole frequency /c. We 
compute fc by measuring the light storage time r = l/(47r/c) in each cavity. 
A single measurement of the storage time is performed by aligning a 



single arm of the interferometer (as in the right panel of Figure A. 17) and 
holding the cavity on resonance using the coil actuators. Then, the power 
transmitted through that arm is recorded as we rapidly take the cavity out 
of resonance. We fit the resulting time series to a simple exponential decay, 
whose time constant is the light storage time in the cavity. This measurement 
is performed several times per arm, and the average light storage time is used 
to calculate the cavity pole frequency. Table [4] shows the values of fc used in 
each model. 
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Table 4: Summary of cavity pole frequencies fc used in each interferometer's sensing 
function in S5. HI and H2 have used the average of each arm, hence their numbers 
reported below are the same, with uncertainty estimated as the quadrature sum of each 
result. 





fc (Hz) 


(Hz) 


HI 
H2 
LI 


85.6 ± L5 
158.5 ± 2.0 
85.1 ± 0.8 


85.6 ± 1.5 
158.5 ± 2.0 
82.3 ± 0.5 Hz 



3.2.3. Analog to Digital Conversion^ ADC{f) 

Each of the four photodiodes used to measure the power at the dark 
port are sampled at 16384 Hz. The dominant frequency dependence of this 
analog-to-digital conversion process arises from the analog anti-aliasing fil- 
ters. These filters are analog eighth order elliptic filters, which differ only 
in corner frequency at the two sites: 7.5 kHz for the Hanford and 8.1 kHz 
for Livingston. The frequency dependence is unity below 1 kHz. Above a 
few kHz, the magnitude changes less than 2%, but the phase loss from these 
filters becomes non-negligible (> 180 deg). The residual frequency depen- 
dence between this model and measured transfer function of the filter is also 
included. The discrepancy occurs only above 1 kHz and varies less than 2% 
in magnitude and 5 degrees in phase. 

3.2.4' Time Dependence^ j{t) 

We measure the time dependence of the sensing function by digitally 

injecting a signal, Sci{f)^ at the output of the digital filters, D{f)^ prior to 
the control signal, s^if)^ at three line frequencies fd near 50, 400, and 1100 
Hz. The time-dependent coefficient ^{t) is defined as 

7(t) = a(t)m = - r ] f ^'''^^i ^ (22) 

where GlUci) is the modeled DARM open loop transfer function at the ref- 
erence time in each epoch at a given calibration line frequency, fd] Sdifd) 
and suifci) are the excitation signal and the control signal, respectively, each 
digitally demodulated at the same frequency and averaged over 60 seconds. 
The coefficient generated from fd ^ 400 Hz is used to scale the response 
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function model; the other two frequencies are used to confirm that the varia- 
tions are independent of frequency. In the ideal case (no noise on top of the 
injected line and with a perfect model for GlUci))-, the coefficient is a real 
factor near unity. Figure 11 shows the evolution of 3?e{7(t)} over the course 
of the science run for each detector. 




Figure 11: Time-dependent corrections to the sensing function 5Re{7(t)} over the course 
of the science run. 

We also separate the relative uncertainty of the time dependent coefficient 
into those of systematic and statistical origin. As the coefficient is 
ideally real and unity, we expect the imaginary part of the measurement 



defined in Eq. (22) to be a random time series with zero mean. A non-zero 



mean would indicate a systematic error in our estimate of 3?{7(t)}, given by 



1 — \^m{'^{t)}j (assuming the real part of 7(t) is unity). The measured 

mean is less than 5% for all detectors, implying a negligible systematic error 
of 0.1% and is ignored. 

The statistical error is determined by the signal-to-noise ratio of the cal- 
ibration line at frequency fcu and is estimated by the standard deviation of 
Q'm{7(t)}, measured in every epoch at a sampling rate of 1 Hz. Though the 
statistical error is roughly equivalent in all epochs for a given detector, we 
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chose the largest standard deviation as a representative error for the entire 
run. Figure 12 shows an example histogram of ^m{'y{t)} for H2. 



Y(t) Uncertainty 
H2, Epoch 1 




lm(Y) 

Figure 12: Histogram of '^m{j{t)} for the epoch 1 in H2, the standard deviation of which 
represents the error estimation for this epoch. 



3.3. Open Loop Transfer Function 

The open loop transfer function, GL{f)^ is measured while the interfer- 
ometer is controlled, operating in the nominal configuration, at designed 
sensitivity. We use a digital DARM excitation with amplitude much larger 
than ALf^xt^ such that we may assume it to be a contribution to measurement 
noise. During the measurement we assume no time-dependent variations oc- 
cur, and set ^{t) = 1. We compare this measurement against our model 
of the open loop transfer function which is the product of each subsystem 
described above (see Eq. [5]), and scale the model by the measurement's mag- 



nitude at the expected unity gain frequency to form ICc as described in ^3.2 



Values for the sensing scaling coefficient averaged over epochs, are shown in 
Table [5l 

We measure the open loop transfer function many times during the course 
of the science run. To compare these measurements against the model for 
each epoch, they are normalized by the magnitude of the open loop transfer 
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Table 5: Average value for scaling coefficients Kc for the sensing function, Cl{^^\) for each 
interferometer. They are stated without uncertainty, since these quantities are derived 
from measurements of the open loop gain and actuation scaling coefficient. See further 
discussion in Q 





Kc {cts/W-^^m) 


HI 


0.15 


H2 


0.61 


LI 


9.1 



function at a fixed unity gain frequency. This normalization removes the 
time dependent scale factors between measurement times such that a fair 



comparison can be made. Figure 13 shows the results of this comparison. 

The uncertainty estimation in the open loop transfer function magni- 
tude and phase {{(J\Gl\/\^l\Y and ) are separated into systematic and 
statistical uncertainty. We expect the ratio of the model and our measure- 
ments to follow a Gaussian distribution with unity mean in magnitude and 



zero mean in phase. This ratio is shown in Figure [T3| We observe a non- 
Gaussian systematic in all detectors from an unknown source, most appar- 
ent in the Hanford detectors. We estimate this systematic uncertainty in 
magnitude and phase by subtracting a smoothed version of the residuals, 
G^^if) = {(J^^^^^ / G^^^^) ^ from unity and zero, respectively. The statistical 
uncertainty, (Jt.\Gl\ crj]^^^^, is estimated from the standard deviation of 
the remaining scatter in the ratio after the systematic error G^£^{f) is sub- 
tracted. Both the systematic and statistical errors are added in quadrature 
to form the total uncertainty in the open loop transfer function model, 

(^)' = (^E|G.|)' + (l-|Gr(/)l)' (23) 

2 



^Icr, = r^^cj (24) 



4. Uncertainty Estimation 

The measurement uncertainty of each component of the response function 
described in ^ are folded into a complex function of frequency known as the 
"error budget." 



28 



S5 H1 DARM Open Loop Gain S5 H2 DARM Open Loop Gain 




Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) 




Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) 



S5 L1 DARM Open Loop Gain 




Frequency (Hz) Frequency (Hz) 




Frequency (Hz) Frequency (Hz) 



Figure 13: Open loop transfer function model vs. measurement comparisons for HI (top 
left), H2 (top right), and LI (bottom) in all of S5. The four panels shown are the magnitude 
and phase of model and measurements (top and bottom left), and the ratio between model 
and measurements (top and bottom right). 



We do not assign any uncertainty to the digital filters D{f) nor directly 
to the time-independent component of the sensing function CL{f). The dig- 
ital filters, which are well-known digital functions, are placed into the model 
without uncertainty. As described in §3.2| the frequency dependence of the 
sensing function is composed of parameters measured to negligible uncer- 
tainty. Uncertainties in its scaling coefficient /Cc, are accounted for in the 
open loop transfer function and actuation function uncertainty. 
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The uncertainties of the remaining quantities in the response function 
A{f)^GL{f)^ and j{t) are treated as uncorrelated. If the uncertainties are 
completely correlated (i.e. there are none in CL{f))^ the covariant terms in 
the estimation reduce the overall estimate of the response function uncer- 
tainty [2^ . Since we do not have an independent estimate of the uncertainty 
in the sensing function, we adopt this conservative estimate. 

We re- write the response function in terms of the measured quantities to 
which we assign uncertainty, 

and separate into magnitude and phase (dropping terms which include the 
uncertainty in D{f)), 



\Rl\ = 



\A\ 



L\ 



l + (7|Gi|)2 + 27|Gi|cos(</,Gj . (26) 



= a-ctan( ^l^f'"i^-' + ^"'';;-^/-' ), (27) 

\^\Gl\ cos {(pA) + COS {(pA - (PGl) ) 

such that the relative uncertainty in magnitude and absolute uncertainty in 
phase are 



G 



\Rj\ 



2 



+ ^m{W^} VJ^^ + ne{Wy ( ^ ) (28) 



2 

7 



L 



+ "SleiWY al^^ + ^m{WY (^^J , (29) 

where we define W = 1/(1 + Gl) [23]. Each uncertainty component in Eqs. 
[28l and [29] is assumed to be the same over the course of the science run 
(independent of epochs). However, the complex coefficient W is different for 
each epoch. 
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Our calculation of the response function includes the open loop transfer 
function model which is approximated by replacing the complete cavity re- 
sponse Hpp{f) (Eq. [9]) with the single pole transfer function Hsp{f) (Eq. 



10) in the sensing function subsystem. We include the ratio of the response 
function calculated with and without the correct cavity response in our error 
budget, 

Rf^if) l + G,if) ^^^^ 

added linearly (as opposed to in quadrature) because the approximation re- 
sults in a frequency-dependent scaling of the response function with known 
sign. As with the weighting function this term involves the direct multi- 
plication of the open loop transfer function and therefore is epoch dependent. 
As an example, Figure [M] shows this error contribution from the third epoch 
in each detector. 



S5 Response Function Error From Single Pole Approximation 




Frequency (Hz) 



Figure 14: Systematic uncertainty in response function arising from the single pole approx- 
imation of the Fabry- Perot cavity response in open loop transfer function. This uncertainty 
is epoch-dependent; only the third epoch for each detector is shown. 
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5. Results 



In Figure [15] we plot the final response function for all interferometers 
for the entire fifth science run. Figure [16] shows the frequency dependence 
of all terms in the error budget of the response function for the third epoch 
of each detector. In Table [6| we summarize the frequency-dependent uncer- 
tainty of each interferometer's response function by dividing the error into 
three frequency bands: 40-2000 Hz, 2000-4000 Hz and 4000-6000 Hz and 
computing the RMS errors across each band, averaged over all epochs. All 
epoch uncertainties are within 1% of the mean uncertainty stated. 



Table 6: Summary of band-limited response function errors for the S5 science run. 





RLif) Mag 
40-2000 Hz 


nitude Ei 

2-4 kHz 


cror (%) 

4-6 kHz 


Riif) Phj 
40-2000 Hz 


ise Error 

2-4 kHz 


(Deg) 

4-6 kHz 


HI 
H2 
LI 


10.4 
10.1 
14.4 


15.4 
11.2 
13.9 


24.2 
16.3 
13.8 


4.5 
3.0 
4.2 


4.9 
1.8 
3.6 


5.8 
2.0 
3.3 



The largest source of systematic error in most data analysis techniques 
used to analyze S5 LIGO data is the uncertainty in response function mag- 
nitude [U El E] • Our inability to measure the sensing function independently 
of the closed loop (specifically its scaling coefficient) forces a conservative, 
uncorrelated treatment of the uncertainty in the measured subsystems, A{f) 
and GL{f)^ inffating the total uncertainty in the response function. In all 
detectors, we find the uncertainty in the actuation function, dominates 
the response function error budget in magnitude. 

The statistical uncertainty in the free-swinging Michelson measurements 
of the actuation scaling coefficient are the primary source of the actuation 
uncertainty. In the Hanford detectors, the uncertainty arises from our in- 
ability to displace the test mass above residual external noise sources at high 
frequency. This decreases the signal-to-noise of the measurement, inffating 
the uncertainty estimate across the measurement band. For LI, in which 
we have obtained a large number of measurements using several methods of 
the free-swinging Michelson technique (see Appendix A), we have found the 
results to be inconsistent with a Gaussian distribution. We attribute this to 
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a poorly understood underlying variation in the technique, for example the 
assumption that the optical gain is time-independent over the course of the 
measurement suite. 

The assumption that the actuation scaling coefficient is linear in ampli- 
tude over the range of actuation, from the 10~^ m employed for the free- 
swinging Michelson technique to the 10~^^ m required to compensate for 
expected gravitational wave signals, has not been confirmed. To investigate 
the linearity of the actuation scaling coefficients over this range of actuation 
amplitudes, and to bound potential overall systematic errors, we have em- 
ployed two additional, fundamentally different, actuator calibration methods. 
The so-called "frequency modulation" technique |19j uses an independently 
calibrated oscillator to frequency-modulate the interferometer's laser light, 
creating an effective length modulation on the order of 10~^^ m while op- 
erating in a single-arm interferometer configuration. The so-called "photon 
calibrator" technique ^] uses auxiliary, power-modulated lasers to displace 
the test masses by approximately 10~^^ m via radiation pressure with the 
interferometer in its nominal configuration (see Figure 1). Both methods are 
employed at select frequencies across the LIGO measurement band. Statisti- 
cal uncertainties for both methods are reduced to the 1% level by averaging 
many measurements. 

At the end of the S5 science run, a detailed comparison between these two 
methods and the free-swinging Michelson technique was performed. With 
all three calibration methods, actuation coefficients were measured over the 
frequency band from 90 Hz to 1 kHz for each end test mass. For the HI and 
H2 interferometers, all calculated actuation coefficients-for all frequencies, 
for all four masses, and for all three methods-were within a ±15% range. 
The maximum difference between the mean value for any method and the 
mean value for all three methods, for any of the four end test masses, was 
3.7% [18j. This indicates that the overall systematic uncertainties in the 
actuation functions determined using the free-swinging Michelson method, 
and therefore the magnitudes of the interferometer response functions, are 
within these bounds. 
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Figure 15: Frequency dependent response function, RlH)^ for the three LIGO interfer- 
ometers for ah epochs of the S5 science run. 
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Figure 16: Frequency dependent error estimates for the response function for HI (top 
left), H2 (top right), and LI (bottom left). In magnitude (top panel) and phase (bottom 
panel), the total uncertainty (dashed-black) is composed of the uncertainty in actuation, 
<^|A|/|^|7 ^4>A (i^ blue), the open loop transfer function magnitude g\q\I\G\ (in green), the 
open loop transfer function phase, (in red), the time-dependent factor, (in cyan), 

and from the single pole approximation of the Fabry-Perot cavity response (in magenta). 
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6. Summary 

The LIGO interferometers have provided the some of the world's most 
sensitive gravitational wave strain measurements during their fifth science 
run. We have described a model used for each interferometer's differential 
arm length control loop known as the length response function, i?L(/, t)^ the 
proportionality between the digital Pound-Drever-Hall error signal and dif- 
ferential displacement of the end test masses. Measurements presented here 
have shown the frequency-dependent uncertainty in i?L(/, t) is less than 15% 
in magnitude and 5° in phase in the frequency band where the interferometer 
is most sensitive. Because we cannot measure the sensing function without 
the interferometers under control, this estimate is limited by our ability to 
measure the actuation function. The results of two fundamentally different, 
high-precision methods for measuring the actuation functions [IHllTTj confirm 
that the free-swinging Michelson results are within the stated uncertainties 

m 

In the two calendar year science run, as our knowledge of the long-term 
characteristics of the instrument increased, a great deal of improvements 
were made to our measurement techniques compared with prior results [7j. 
However, future detectors will have more sophisticated actuation and sensing 
methods [20]. In addition, an amplitude uncertainty of 10% or less is 
required to reduce the calibration uncertainty below other systematic errors 
in the continually improving astrophysical searches [26j. To achieve this 
goal, the non-Gaussian distribution of the actuation function measurements 
must be better understood and independent techniques of measuring the 
actuation coefficient, like the frequency modulation and photon calibrator, 
must be used in concert with the standard techniques presented in this paper 
to reduce limiting systematic errors. 
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Appendix A. The Free-Swinging Michelson Techniques 

The technique used for determining the actuation coefficients, /C^'^, for 
the fifth science run is known as the "free-swinging Michelson" technique. 
This technique uses the interferometer's well-known Nd:YaG laser wavelength 
(A = 1064.1 ±0.1 nm, [24^ 36j) as the calibrated length reference while 
using the test mass coil actuators to cause a change in length of simple 
interferometer configurations. The technique may be used in two similar 
methods: the "Simple Michelson" and "Asymmetric Michelson" methods. 

The Simple Michelson method is composed of two steps. The first step de- 
termines the actuation scaling coefficient for the input test masses Ki with the 
interferometer in a non-standard configuration called a frequency-modulated 



simple Michelson (see left panel of Figure A. 17). The second step determines 
the end test mass actuation coefficient, /C^, from the input test coefficient, 
/C^, and transfer function measurements of the input and end test masses of 



a single Fabry Perot arm cavity (see right panel of Figure A. 17). The Asym- 



metric Michelson determines ICa directly using the configuration shown in 



Figure A. 20. Both free-swinging Michelson methods are described below. 



Simple Michelson Method 

The Pound-Drever Hall error signal at the anti-symmetric port, qAS) for 
a frequency-modulated simple Michelson interferometer is 

1 /47r \ 

qAs = -App sin f y A^j , (A.l) 
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Figure A. 17: Interferometer configurations used during the simple Michelson method of 
measuring the actuation scaling coefficient. Left: The simple Michelson configuration, 
where the power recycling mirror and end test masses are misaligned. Right: The single 
arm configuration, with the power recycling mirror and the opposing arm's input and end 
test mass are misaligned. 

where App is the peak-to-peak amplitude of the signal (proportional to the 
input power, the product of Bessel functions of modulation strength, and the 
transmission of the sidebands into the antisymmetric port from the Michelson 
asymmetry), A = 1064 nm is the wavelength of the input laser light, and 
A£ = ix ~ is the differential arm length of the Michelson. App is measured 
by aligning the simple Michelson and recording the qas time series as it is 
left uncontrolled. In this configuration, external noise sources (e.g. residual 
ground motion) are large enough to cause the Michelson to sweep through 
many interference fringes. 

For the simple Michelson, when A£/X <C 1, 

qAs - kA£ (A.2) 

with the simple Michelson's "optical gain," 

k = i27r/X)A„, (A.3) 

which has units of digital signal counts per meters of input test mass motion. 
After a measurement of App is obtained, we control the optics using their 
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coil actuators, forcing the Michelson into the hnear regime where Eq. |A.3| is 
vahd. 

The actuation function of the suspended input test masses can be ap- 
proximated by the center-of-mass force-to-displacement transfer function, 
Pl^ with a scaling coefficient, /C^. We obtain a measurement of Ki for a 
given input test mass by introducing a digital excitation exci into the control 
loop that is much larger than residual external noise sources. The excitation 
is performed over many frequencies in the gravitational wave band; assum- 
ing the model is complete, the coefficient should be frequency-independent 
across the band. We obtain a solution for the digital excitation counts on 
the input test mass in terms of meters of resulting motion as measured by 
qas (normalized by the pendulum response Pcm)^ 

- - m (^) (i) ^ 

The first term is the measured response of the Michelson during the single 
input test mass excitation. The second term contains the open loop transfer 
function Gsm of the simple Michelson control loop (measured just prior to 
measuring the response to excitation) and the quantity k is as defined in 



Eq. A. 3, We take the median of Ki (denoted with "bra" "kets," ( )), over 
the measured frequency points to remove measurement outliers and residual 
frequency dependence (or time dependence of /c, as discussed in ©1). Figure 



|A.18| shows an example measurement of Ki for each input test mass in H2. 



We then configure the interferometer to form a single Fabry-Perot cavity 
composed of one arm of the interferometer, and control it such that the cavity 



is under resonance (see Figure A. 17). In this configuration, the response of 
the single arm cavity (now recorded by the in-phase demodulated output 
lAS) see [35] for details) to sequential length excitations of the input test 
mass, exQ, and end test mass, exCg, are measured. The ratio, IZie of these 
to transfer functions can then be used to write the actuation coefficient for 
the end test masses as 

. ^„ (S.) . (i^) (££^) (S.) 

where Ka has units of test mass motion in meters (as measured by qAs) per 
count of digital excitation. Figure [A. 19| shows a measurement of Ka for each 



test mass in HI. As in the first step, the median of the frequency points 
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Figure A. 18: Example actuation scaling coefficients for the H2 input test masses /Cf'^, 
measured using the simple Michelson method. Top: JCi as a function of frequency for the 
X (red) and Y (blue) input test mass. Solid lines indicate the median of the data points 
{JCi), dashed lines indicate la error bars (the standard deviation of all frequency points). 



measured in /C^ is used to form a single value for the coefficient over the 
measurement bands. 

Asymmetric Michelson Method 

During the latter part of the science run, a more direct approach of deter- 
mining the actuation coefficient JCa was taken, using the "asymmetric Michel- 
son technique." This method is similar in principle to the simple Michelson 
version of the free-swinging Michelson technique, however, we configure the 
interferometer as shown in Figure |A.20[ In this method the response of the 
end test mass is measured directly and 



k \R 



cm 



Figure [A.21| shows an example result for LI. 

The quantity k may vary slowly over the measurement period due to in- 
put laser power fluctuations, interferometer alignment, etc. The asymmetric 
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Figure A. 19: Example actuation scaling coefficients for the HI end test masses /C^'^(/), 
measured using the simple Michelson method. Top: JCa as a function of frequency for the 
X (red) and Y (blue) end test mass. Solid lines indicate the median of the data points 
(/C^'^), dashed lines indicate la error bars (see §3 for description). 



Michelson is particularly sensitive to these variations as round trip power 
loss is large. For this method, we employ a more sophisticated technique for 
determining the amplitude App^ developed originally by RoUand et al [32j. 
A plot of qas versus the total power incident on the photodiodes should be 
an ellipse whose semi-minor axis is App/2. We obtain a fit to this ellipse and 
extract App with a quantifiable statistical error. 

Results 

Using the above methods, the actuation coefficient is measured many 
times for each optic in each interferometer over the course of the science 
run, and their mean used as the actuation scaling coefficient for all model 
epochs. Figure |A.23| shows the representative median and estimated uncer- 
tainty for each of these measurements. Table A. 7 summarizes the actuation 
coefficients used in the actuation model, (JCa) for the three interferometers in 
the fifth science run, using either simple Michelson or asymmetric Michelson 



techniques, with statistical uncertainty as described in ^3.1.4 
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Figure A. 20: "Asymmetric Michelson" configuration of the interferometer. With the power 
recycling mirror misaligned, an input test mass and opposing end test mass are aligned. 

Table A. 7: Summary of the actuation scaling coefficients measured during S5. These single 
numbers are formed by the mean of each measurement's median {JCA)j (6 for each end 
test mass in HI, 5 in H2, and 14 and 15 for the X and Y test masses, respectively in LI). 
Only statistical uncertainty is reported here; systematic uncertainty is folded the the total 
uncertainty of the actuation function. 





/C^ (nm/ct) 


/C^ (nm/ct) 


HI 
H2 
LI 


0.847 ±0.024 
0.934 ± 0.022 
0.434 ± 0.039 


0.871 ±0.019 
0.958 ± 0.034 
0.415 ±0.034 



Each simple Michelson measurement of a given optic's coefficient is as- 
signed magnitude and phase uncertainty, 



)■ = 


"std(|/C,|)" 


2 


1 Std(|7^^e|)1' 




[ m) \ 


± 


.(|7^^e|) ^/N 



a 



Std(0^j2^ 



std(07e,. 



(A.7) 
(A.8) 



and asymmetric Michelson measurement is assigned magnitude and phase 
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Figure A. 21: Example actuation scaling coefficient for the LI X arm end test mass /C^, 
measured using the asymmetric Michelson method. Top: ]C\ as a function of frequency. 
Solid lines indicate the median of the data points (/C^), dashed lines indicate la error 
bars. Bottom: Histograms for /C^, with the median shown in solid pink, and error bars in 
dashed pink. 



uncertainty 



a 



std(|/C^|) 
. {\}Ca\) 



(A.9) 
(A.IO) 



In simple Michelson technique, measurements of JCi were found to be 
inconsistent with a Gaussian distribution across the frequency band. We 
therefore estimate the uncertainty in the median, (/C^) to be the standard 
deviation alone. However, in the second step (Eq. |A.5[ ), we have found the 
single arm transfer function ratio, i?^e, to be consistent with a Gaussian dis- 
tribution across the frequency band, so we estimate the median uncertainty 
as though it were a gaussian distribution and divide the standard deviation 
by VN where N is the number of frequency points. In the asymmetric Michel- 
son method, where the measurement of JCa is similar to that of ICi in the 
simple Michelson method, we again do not assume a Gaussian distribution 
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Demodulated Signal vs. Total Power at Antisymetric Port 
Free-Swinging Asymmetric Michelson; L1 Oct 27 2007 




Figure A. 22: An example ellipse produced by photodiode demodulated signal qas versus 
total power, and corresponding fit to the ellipse used to determine App in asymmetric 
Michelson method. 



over the measurement band, and take the standard deviation alone. 
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Figure A. 23: Measurements of the actuation scaling coefficient /C^'^, measured over the 
course of the fifth science run. Measurement numbers 6 in HI and 7 through 15 in LI 
used the asymmetric Michelson technique, the remainder were measured with the simple 
Michelson technique. (Left) Magnitude and phase of each measurement median {JC\)j 
(Top) and (/C^)j (Bottom) treated independently. The statistical uncertainty of actuation 
function is the quadrature sum of each arm's uncertainty, which takes the larger of the 
standard deviation of each measurements median, {JCA)j or the mean uncertainty divided 
by the number of measurements (JKa^j- 
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